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SUMMARY 



A procedure for the shear— lag analysis of box teams, 
such as wing structures, is outlined. A previously pub- 
lished method separated the most essential part of the 
shear-lag analysis from the Kc/I analysis. The present 
method, "by slightly modifying the computation of the chord- 
wise distribution of stress, entirely separates the shear- 
lag analysis from the He/ I analysis. The discussion points 
out that savings in time and greater accuracy will result 
from this separation. 



A comprehensive paper on shear lag (reference 1) pub- 
lished previously gives two methods of analysis. The 
method recommended for ordinary use achieves some of its 
simplicity and sueed by separating the main part of the 
shear-lag analysis from the Mc/ I analysis; the two anal- 
yses are combined just before the chordwise distribution 
of the stringer stresses is computed. When the comuutation 
of the individual stringer stresses is considered necessary, 
the method recommended in reference 1 should be slightly 
modified to keep the shear-lag analysis separate from the 
He/ I analysis until the final addition of stresses is made. 
The present paper gives the complete outline of the modi- 
fied procedure and discusses the advantages gained by sep- 
arating the shear-lag analysis from the Mc/l analysis. 
The treatment is chiefly intended for wing structures of 
the central— box type. 



INTRODUCTION 



♦ 



THE SHEAR-LAG ANALYSIS OP BOX BEAMS 



General Procedure of Analysis 
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The general procedure for shear-lag analysis consists of 
the following stepsj 



1. At a number M,.jtat ions along the span, prelim- 
inary section modulilkre calcuriat ed ; estimated amounts of 
skin are assumed to work with the stringers on the com- 
pression side. Preliminary values of the stringer stresses 
are calculated with these moduli. 

2. The preliminary stringer stresses are used to cal- 
culate the effective widths of skin. Marguerre 's< formula 
for the effective width 

2w/b = Jo CT lo (1) 

is recommended as a good compromise between accuracy and 
simplicity. Average values of the critical stress a cr 

of the sheet and of the stringer stress a may he used 
for any station unless the thickness of the skin varies 
chordwise .. 

3. After the effective widths of sheet have been 
established by checking them for agreement with the as- 
sumed values, final section moduli are calculated and the 
final He/ I analysis is made. 

4. The shear— lag analysis is made by the methods to 
be described later; the result of this analysis is a set 
of stress corrections. 

* 

5. The stress corrections resulting from the shear- 
lag analysis are added to the stresses resulting from the 
final Me/ I analysis. 

It should be noted that .steps 1 to 3 constitute the 
He/ 1 analysis familiar to every analyst, with the excep- 
tion, perhaps, of formula (1) for the effective width. 
There are no restrictions placed on the Mc/l analysis; 
lift and drag forces may be considered separately or 
jointly and principal axes may or may not be used, as de- 
sired. 



Basic Procedure for Shear-Lag Analysis 

The basic procedure for shear-lag analysis is the 
procedure applicable to box beams of rectangular cross 
section symmetrical about a spanwise plane parallel to the 
plane of the loading being analyzed. As a rule, only beam 
loads need be investigated; the shear-lag effect on chord 
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loads may be neglected. Lift loads may often "be used in 
place of team loads. Symmetry being assumed to exist, all 
calculations are made for the half-section (fig. 1). 

q The individual steps of the basic procedure are as 

-■j- follows; 

1. The beam is divided into a number of bays, so that 
the properties of the cross section and the applied run- 
ning shear may be considered as sensibly constant within 
each bay. The stations are numbered as shown in figure 2. 

2. For each bay, the following properties are com- 
puted : 

Ay area of corner flange- (This area consists chiefly 
, of two items: the area of the corner angles when 
they exist, and an area l/6 ht^ that takes into 
account the portion of the bending carried by the 
shear webs . ) 

A L area of all longitudinals (stringers and effective 
widths of skin) 

A T = Ay + A L 

b. width of half-section 
b r width at root 

t thickness of cover sheet (chordwise average if there 
is a variation) 

tf fictitious thickness of cover sheet (t b r /b) 

0 effective shear modulus of cover sheet (This modulus 
should be estimated on the basis of the combined 
stresses acting on the sheet at the load being in- 
vestigated — for instance, limit load or ultimate 
des ign load. ) 



K 



shear-lag parameter j K s = "eTT ^ r + A/ J ^ 



L length of bay 
Gt f tanh KL 



p 3 3L_ (3) 



4 



_K_ 

q Gt f sinh KL 



(4) 



V = (5) 
Gtf h A j 

where is the shear force in the shear web; that is, 

the external shear force S ? reduced by the vertical com— 

ponent of the stringer forces if the beam tapers in depth 

S W = H - I 8 (6) 

For numerical work, it will "be found convenient to 
modify the expressions f or p , q , and Y "by writing G 
in terms of a reference modulus — for example, the shear 
modulus of the material. The value of the reference modu- 
lus may then "be omitted from the expressions for p, q, 
and Y, as will he seen from an inspection of the next 
step. The thickness t may he treated in the same manner. 

3. A set of equations is written for the statically 
indeterminate forces X, in which the coefficients p, q, 
and Y are used, 

* 0 <li - XjCp! + p 3 ) + *2q 2 = - Y x + Y 2 
X^s - X 2 (p 3 + p 3 ) + X 3 q 3 = - Y 2 + Y 3 



(7) 

'n-i^n " X n<Pn + Pn+ i> + X n+ i*n+ i * ~ V n + Y n+! 



X r -iq r - x r(Pr + Pr+i) ■ ~ Y r + Y r+i 

4. Some of the quantities appearing in the set of 
equations are determined hy the boundary conditions as 
f ollows : 

(a) If station zero is at the extreme tip, X 0 = 0. 

(b) If station zero is at the joint between the end 

of the box beam and a special tip structure 
connected in such a way that only the corner 
flanges carry stress across the joint, 



(c) If the longitudinals stop at the root section 
and do not carry stress into the fuselage, 
the last equation of the system drops out 
and 

X„ = ^ (9) 
l T 



r hAr 



In formulas (8) and (9), the terms on the right- 
hand side are , of course, those ap-propriate 
to the section under consideration. 

(d) If the longitudinals do not stop at .the root 
section; that is, if the wing cover being 
analyzed is continuous through the fuselage, 

r+i v 

because S„ = 0 in the fuselage and 

= + < l0 > 

5. The set cf equations is solved and the values of 
X are tabulated. 

6. The correction to the Mc/ I stress in the corner 
flange at a given station is 

Aer F = (11) 
Ay 

The correction to the He/ 1 stress, in a stringer located 
at a distance y from the center line of the beam is 

La = A Op - D [l J (12 > 



where 



D = t Tl(±- + IS) (13) 
3 \k F k L J 



The addition of the corrections given by formulas (11) and 
(12) to the Mc/I stresses completes the analysis of the 
stringer stresses (fig. 3). The shear stresses, in the 
sheet may then "be obtained by statics from the equilibrium 
of longitudinal forces on sections of the cover. 

For a single concentrated load applied to a beam of 
constant cross section, figure 4 shows the Mc/ I stresses 
in the flange, the shear-lag corrections AOf = X/ Aj> , and 

the total flange stresses. It will be noted that the sign 
of the correction changes along the span and that there is 
one station at which there is no correction to the Mc/ I 
stress. Similar curves apply to beams with distributed 
loading. The station with zero correction is of practical 
interest when a strain survey is to be made at a single 
station, because strain readings taken too close to this 
station would be misleading. 

The formula (12) is an approximation ; that is more 
convenient for practical use than the theoretical expres- 
sion involving hyperbolic functions (reference l), The ap- 
proximation breaks down to some extent near stations where 
the longitudinals stop. In such regions, the stringer 
stresses near the center line approach zero, while the 
stresses given by the approximation (12) pass through zero 
and change sign. No practical significance attaches to 
this failure cf the approximation, because it occurs only 
where the stresses are too small to be of interest. 



Modifications of Basic Procedure 

Utilisation s_f or _camber_of _c over^— Only slight modi- 
fications are necessary to adapt the procedure to the 
analysis of beams with slightly cambered covers. The ex- 
pression (2) is replaced by 



where b' is the developed half— width of the cover sheet 
and c x and c s are the effective cambers, that is, the 
vertical distances from the corner flanges to the centroids 
of the stringer forces. The stresses being unknown at the 
outset, the locations of the force centroids must be esti- 
mated; for most practical purposes, the centroids of the 




(14) 
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stringer areas Aj, may be u*ed (fig. 5). The formulas 
(5) and (6) are replaced, by the more general expressions 

V * $ Jtt^k (15) 
Gt f I 

* S E - tan S (16) 

where a is the distance from the neutral axis and Q is 
the static moment of A^ and Aj. about the neutral axis. 

From this poi-nt on, the procedure fcr calculating the; X- 
forces is identical with the basic procedure. 

The correction to the Mc/ 1 stress of the corner 
f lange is 



where c is the effective camber of the- cover being anal- 
yzed. The correction to the Mc/l stress in a stringer 
located at a distance y from the center line of the beam 
is given by the expression 

:: '.: "r f 1 -©"] (18) 

where D is given by formula (13). 

The shear-lag effect decreases the stringer stresses 
near the center line of the beam and increases them near 
the corner flange. This shift decreases the effective 
beam depth when the. cover is cambered; a correction must 
therefore be applied not only to the side being analyzed - 
for example, the compression side - but also to the oppo- 
site side, the tension side. The magnitude of the correc- 
tion force is Xc/hy. This force increases the flange 
stress on the opposite side when X is positive. 

Modification^ for _un.sy_mme_tr leal _er qss sections .— No 
satisf actory~theoretical treatment of shear-lag effect in 
beams with unsymmetr ical cross sections appears to have 
been published so far, There is in the literature a paper 
dealing with the analysis of beams of unsymmetr ical closed 
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section and a paper dealing with beams of unsymmetr ic al 
open cross section. Both papers contain several assump- 
tions of questionable validity. The second paper, more- 
over, omits the vertical component of the shear in the 
cover sheet; at best, then, the formulas given in this 
paper would be approximations valid only when the dissym- 
metry is very small. 

The lack of a satisfactory theory makes.it impossible 
at present to give a rational procedure for beams with 
very unsymmetr ical cross sections. Fortunately, the dis- 
symmetry of practical cross sections is seldom very marked, 
and it is permissible to resort to the common expedient of 
using a mean cross section obtained by averaging the val- 
ues of h w and Ay for the front and rear spar6. It 
should be carefully noted that th e mean cross section is 
used only to calculate th,« ajM> correct I ens , Bot__ th* 
M c / I stresses. The error committed by using a mean sec- 
tion then affects only the shear-lag corrections ha and 
is, therefore, a very small fraction of the total stress, 
the shear— lag correction itself being only a fraction of 
the t otal stress. 

Aaalyai.s_£l_li£.affifi.--Kii.li_c.u.i = 5iii£. Jt - The analysis of 
beams with small or medium cut-*outs is effected conven- 
iently by the method of liquidating forces (reference 1). 
The beam is first analyzed on the assumption that no cut- 
out has been made. The internal forces at the boundaries 
of the proposed cut-out are. next calculated. Liquidating 
forces are then introd^lced, equal and opposite to the 
boundary forces just found. The stresses caused by the 
liquidating forces are calculated and are superposed on 
the stresses calculated for the box without cut-out. The 
details of the calculation for certain typsB of cut-outs 
may be found in references 2 and 3. 

Beams with full-width, or nearly full-width, cut-outs 
are analyzed in two parts. The part outboard of the cut- 
out is analyzed as a beam with longitudinals interrupted 
at the "root," that is, at the outboard end of the cut-out. 
The part inboard of the cut— out is analyzed as a beam with 
a wing tip structure introducing moments at the tip, that 
is, at the inboard end of the cut— out . 



Preliminary Estimates of Shear-Lag Effects 

The general procedure for analysis is very flexible 
and can be easily adapted to the requirements of preliminary 
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analysis or first estimates. The Me/l analysis, as well 
as the shear— lag analysis, can he simplified "by subst itut — 
ing a mean surface without camber for a cambered surface. 
Tho. shear-lag analysis can he simplified by reducing the 
number of bays into which the beam i8 divided. 

For a beam with longitudinals continuous through the 
fuselage, a first estimate of the shear— lag effect may be 
made under the assumptions thpt the carry— through members 
are rigid, that the load is uniformly distributed along 
the span, and that the cross section is constant. The 
theoretical formulas for this case give for the ratio of 
the actual flange stress at the root to the Mc/l stress 



+ 4~ (19) 

and for the ratio of th* average stringer stress (not the 
stringer stress at the center line) to the He/ I stress 

where L is the semispan. When a concentrated load is 
applied at the tip, the factor 2 appearing in formulas 
(19) and (30) must be omitted. This difference must be 
borne in mind when the formulae are applied to test beams 
loaded at the tip or when empirical coefficients derived 
from such test beams are applied to wings. ' 



DISCUSSION 



The most important feature of the suggested method of 
analysis is the complete separation of the shear— lag anal- 
ysis from the Mc/ 1' analysis. In contrast with the meth- 
ods that do not effect such a separation, the following 
advantages are gained? 

1. The simplifying assumptions that are necessary at 
the present state of knowledge (for instance, the assump- 
tion of symmetry of the cross section) cause errors in 
only the stress corrections Aa , not in the Mc/ I stress^. 
The stress corrections being small compared with the Hc/l 
stresses, the errors will be very small compared with the 
t otal stresses , . 
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2. The time requirpd is reduced "by "breaking the anal- 
ysis down into independent operations (Mc/ I analysis and 
shear-lag analysis). This "breakdown makes it possible to 
employ more men simultaneously, in accordance with a recog- 
nized principle of production. 

3. The "breakdown into independent operations makes it 
possible to use each man more efficiently. The engineer 
capable of dealing with shear-lag analyses is not required 
to make routine He/ I analyses. 

4. The method is very flexible and can be easily: 
adapted to rough estimates, preliminary analyses, or final 
analyses by varying the number of the bays for the shear- 
lag analyses and by using well— known short-cuts for the 
Mc/ I analyses . 

■ Other features of the method that should be taken in- 
to account when making comparisons with other methods are 
the following: 

No empirical coefficients are used. This point is 
important because shear lag is affected by so many vari- 
ables that it is very difficult to deduce reliable empir- 
ical coefficients from tests. Factors that are often com-, 
pletely neglected - for instance, the properties of the 
carry-through members in the fuselage — may decisively 
influence the stresses in a given case. Empirical coef- 
ficients obtained by neglecting such factors may be very 
mi s leading . 

The difficulty of obtaining reliable experimental 
data is accentuated by the fact that structures built up 
from sheet metal show considerable irregularities in 
behavior. Extensive strain surveys on a large number of 
specimens should therefore be made before conclusions are 
drawn about the reliability of a proposed method of anal- 
ysis. A small number of strain measurements on a single 
wing might be acceptable for proving the validity of a 
theory that does not contain any questionable simplifying 
assumption; the more questionable the simplifying assump- 
tions are and the more the method relies on empirical 
coefficients, the more necessary it becomes to extend the 
range of the experimental proofs. On the basis of experi- 
mental proofs submitted, the method given here scores over 
any other published method. 

Only mathematical processes with which the average 
analyst is entirely familiar are employed. The most dif- 
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ficult step is the solution of a system of equations anal- 
ogous to the ordinary three-moment equations. The average 
analyst is acquainted with such equations; he is not well 
acquainted, on the other hand, with Fourier series or with 
the solution of differential equations. 

In conclusion, it may be pointed out that every step 
of the analysis is a calculation performed in a definite 
manner. The analyst is neither expected nor required to 
substitute engineering judgment for calculation. Engineer- 
ing judgment is indispensable in design and can be effec- 
tively used to speed up preliminary analysis. For obvious 
reasons, however, a method of final analysis should be as 
free as possible of personal factors. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for AeronauticB, 
Langley Field, Va. . 
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Figure 5.- Box beam with cambered covers. 



